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Forward Pion Production in Hadron Collisions at STAR 

G. Rakness a , for the STAR Collaboration 

a Pcnn State University/Brookhaven National Laboratory, 
Physics 510-A, Upton, NY 11973 

Measurements are reported of the production of high energy tv° mesons at large pseudorapidity, coincident with 
charged hadrons at mid-rapidity, for proton+proton and deuteron+gold collisions at ^/sjvjv = 200 GeV. The p+p 
cross section for inclusive 7r° production follows expectations from next-to-leading order perturbative QCD. A 
suppression of the back-to-back azimuthal correlations was observed in d+Au, qualitatively consistent with the 
gluon saturation picture. Experimental uncertainties regarding the inclusive measurement are discussed. 



1. INTRODUCTION 

In contrast to the case for the nucleon, lit- 
tle is known about the distribution of gluons 
in nuclei, especially in the region of Bjorken- 
x < 0.01, where x is the fraction of the nucleon's 
momentum carried by the parton. Studies of 
deuteron(proton)+nucleus collisions at large cen- 
ter of mass energy (y / siviv=200 GeV) can provide 
constraints on the gluon density in heavy nuclei. 
In the perturbative QCD explanation of parti- 
cle production at large pseudorapidity, a large- 
x parton scatters from a low- a; parton and then 
fragments into the observed particle(s). Forward 
charged particle production is found to be sup- 
pressed in d+Au collisions pQ, consistent with the 
expectation of gluon saturation |2I3| . possibly in- 
dicating a different mechanism for particle pro- 
duction. Explanations of the suppression based 
on leading-twist pQCD calculations using a model 
of gluon shadowing have also been suggested 0]. 
Further tests of the possible role played by gluon 
saturation at RHIC energies are provided by the 
study of particle correlations |5l(i| . 

At ^/s=200 GeV and larger collision energies, 
there is quantitative agreement between NLO 
pQCD calculations and measured p+p cross sec- 
tions at mid-rapidity 7 . This agreement has 
been found to extend to 7r° production at (77) = 
3.8 Further tests of the underlying dynamics 
responsible for forward particle production can 
be obtained from the study of particle correla- 



tions. In particular, strong azimuthal correlations 
of hadron pairs are expected when particle pro- 
duction arises from 2 — ► 2 parton scattering. 

This paper reports cross sections for for- 
ward inclusive ir° production for p+p colli- 
sions at -^=200 GeV. The azimuthal correla- 
tions between a forward ir° ((rj 7T )=A.0) and mid- 
rapidity charged hadrons were studied. In addi- 
tion, exploratory studies with d+Au collisions at 
y / s]v7V=200 GeV are reported, and azimuthal cor- 
relations of hadron pairs are compared to those 
for p+p collisions. The inclusive yields of ir° 
mesons in p+p and d+Au collisions will be forth- 
coming, and details pertinent to this analysis are 
presented here. Forward tt° production in d+Au 
collisions refers to observation of the ir° in the 
direction of the incident deuteron. 



2. EXPERIMENT 

The Solenoidal Tracker At RHIC (STAR) is 
a multipurpose detector at Brookhaven National 
Laboratory. One of its principal components 
is a time projection chamber used for tracking 
charged particles produced at |f?| <1.2. A for- 
ward 7T° detector (FPD) comprising 7x7 matrices 
of 3.8 x 3.8 x 45 cm 3 Pb-glass detectors (towers) 
was installed ~ 800 cm from the interaction re- 
gion near the beam pipe to detect high energy 
7T° mesons with 3.3 < 77 < 4.1. Data were col- 
lected over two years of RHIC operations. In the 
2002 run, p+p collisions were studied with a pro- 
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Figure 1. Inclusive ir° production cross section 
versus ir° energy (E„) at average pseudorapidi- 
ties ({if/) 3.3 and 3.8. The inner error bars are 
statistical, and are smaller than the symbols for 
most points. The outer errors combine these with 
^-dependent systematic errors. The curves are 
NLO pQCD calculations evaluated at T) = 3.3 and 
3.8 using different fragmentation functions. 



totype FPD [5]. In the 2003 run, p+p collisions 
were studied and exploratory measurements were 
performed with d+Au collisions. 

3. DIFFERENTIAL CROSS SECTION 

The differential cross section for inclusive tt° 
production for 30 < E„ < 55 GeV at (77) = 3.8 
was previously reported [Hj- The event recon- 
struction and normalization methods were ex- 
tended to allow measurement of the differential 
cross section at (77) = 3.3 9 . The results are 
shown in Fig. ^hi comparison to NLO pQCD cal- 
culations evaluated at ?7=3. 3 and 3.8. The NLO 
pQCD calculations are consistent with the data, 
in contrast to 7r° data at lower y/s . The inclu- 



Figure 2. Diphoton invariant mass spectra as 
a function of energy deposited in the calorime- 
ter for d+Au data and HIJING+GEANT simu- 
lations. The histograms are reconstructed simu- 
lation events. The points are data with statistical 
errors, normalized to equal area in each bin. 



sive yield was also found to agree with predictions 
of the PYTHIA event generator JT], which uses 
initial- and final-state parton showers to model 
effects beyond leading order. 

Nuclear effects on particle production are quan- 
tified by RdAu j the ratio of the inclusive yield of 
particles produced in d+Au collisions to p+p col- 
lisions. The value of RdAu for forward ir° mesons 
is expected to be smaller than was seen for nega- 
tive hadrons £Q due to isospin suppression of the 
p + p — > h~ + X process 221- We proceed to 
discuss considerations relevant to an accurate de- 
termination of RdAu for forward ir° production. 

The data in Fig. ^ show that the cross section 
falls rapidly with E v . Accurate yield determi- 
nations require accurate energy calibrations and 
good energy resolution. The energy calibration 
of the FPD is performed by reconstruction of the 
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Figure 3. Simulation study of the fractional 
difference of reconstructed and simulated energy 
versus simulated ir° energy. The energy calibra- 
tion is determined to an accuracy of w 1%. Sys- 
tematics of reconstructing the di-photon opening 
angle are seen beginning around 50 GeV. 



invariant mass of the tt° meson, M 77 . An analy- 
sis of the topology of the energy deposition, based 
on measured shower shapes |13| . is used to recon- 
struct the relative energy of the two photons and 
their separation (cZ 77 ). The ir° energy is taken to 
be the total energy in the calorimeter. The gain 
of each tower is determined iteratively by associ- 
ating M 77 with the tower containing the most en- 
ergy in the event. After convergence, an energy- 
dependent correction is applied to account for 
material between the interaction region and the 
detector, which causes decay photons to shower 
prior to reaching the FPD. It also corrects for the 
non-linearity implied by using digitizers with a 
finite number of bits. The energy can be accu- 
rately associated with M 77 up to a point set by 
the systematics of reconstructing c? 77 . At 50 GeV, 
rf 77 ps 1.1 X the linear dimension of a tower. The 
M 77 distribution can be seen in Fig. with res- 
olution < 20MeV/c 2 from 30 - 55 GeV. 

A simulation is used to model the events us- 
ing the PYTHIA OH JING) [Ti| generator for p+p 
(d+Au) collisions, together with a GEANT sim- 
ulation to model the detector response, including 
intervening material. After calibrating the simu- 
lation with the same technique as the data, the 
overlay of data and simulation for M 77 can be 
seen in Fig. [21 The simulation describes the data 
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Figure 4. Coincidence probability as function 
of azimuthal angle difference between a forward 
7T° and a midrapidity leading charged particle for 
p+p (left) and d+Au (right) collisions. The up- 
per plots are simulation using PYTHIA 6.222 and 
HIJING 1.381 described in the text, while the 
lower plots are data. 



well for M 77 as well as several other kinematic 
variables. The simulations are used to determine 
the efficiency of 7r° detection. As was observed in 
P3, the detection efficiency is dominated by the 
geometrical acceptance of the calorimeter. 

Using the simulation, a comparison of the kine- 
matics of events that contain ir° mesons with re- 
constructed variables gives an estimate of the ac- 
curacy of the analysis techniques used. The frac- 
tional difference in energy between the generated 
and reconstructed ir° meson in bins of generated 
7T° energy can be seen in Fig. |3J The uncertainty 
on the energy calibration is « 1%, implying an 
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uncertainty on the yield of w 10%. 

The data in Fig. 2] show that the cross section 
increases rapidly with r\. Accurate yield determi- 
nations require accurate position determinations. 
A Beam Position Monitor (BPM) is located in 
the vicinity of the FPD, and has been accurately 
surveyed. Relative to the BPM, the absolute po- 
sition of the FPD is presently measured with an 
accuracy of 5 mm, giving an absolute accuracy of 
drj = 0.02 at r\ = 4.00. This implies an uncer- 
tainty on the yield of 11%. 

4. CORRELATIONS 

Correlations between a ir° produced at large 
rapidity with large Feynman x and charged par- 
ticles produced at midrapidity were studied with 
p+p collisions and d+Au collisions. Details of 
the analysis, including detailed comparisons to 
simulations, may be found in Ref. The 
leading charged particle (LCP) analysis selects 
the mid-rapidity track with the highest px >0.5 
GeV/c, and computes the azimuthal angle dif- 
ference A(f> — <p^a — 4>lcp for each event. The 
normalized A<p distributions are fit with the sum 
of a constant and a Gaussian distribution cen- 
tered at Acf> — 7r. The fit parameters are the 
area under the Gaussian (S) , representing the az- 
imuthally correlated coincidence probability; the 
uncorrelated coincidence probability (B); and the 
Gaussian width (as)- The normalized Acj) distri- 
bution for d+Au collisions is shown in comparison 
to p+p for both simulations and data in Fig. 0] 
The simulations account for detector resolution 
and reconstruction efficiency for both the forward 
7T° and the midrapidity charged particles. 

We observe a large increase of BdAu relative 
to B pp . The growth in B arises from additional 
nuclcon-nucleon collisions as the deuteron and 
Au beams interact. For the data, a s Au is much 
smaller than <7g P , most likely reflecting the in- 
adequacy of the functional form used to repre- 
sent A^dAu- The azimuthally correlated 7r° + h 
coincidence probability is smaller in d+Au colli- 
sions than in p+p collisions, qualitatively consis- 
tent with behavior predicted in a gluon saturation 
picture UJ. Complete assessment of systematic 
errors is underway. 



5. SUMMARY 

In summary, cross sections for the inclusive 
production of 7r° mesons in p+p collisions at 
y / s=200 GeV, at (77,r)=3.3 and 3.8 are consistent 
with NLO pQCD calculations. The azimuthal 
correlation between pairs of hadrons separated by 
large Arj is described by PYTHIA , which uses 
leading order pQCD with parton showers. Agree- 
ment of these calculations with the inclusive cross 
section and di-hadron correlations suggest that 
forward ir° production arises from partonic scat- 
tering at this collision energy. Exploratory stud- 
ies of forward 7r° production in d+Au collisions 
suggest that the azimuthally correlated compo- 
nent of hadron pairs separated by large Ar) is sup- 
pressed relative to p+p collisions. More data for 
forward particle production and di-hadron corre- 
lations in d+Au collisions are required to reach a 
definitive conclusion about the possible existence 
of gluon saturation in the Au nucleus. 
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